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oo '. Abstract 
O 



Recent advances in space-qualified atomic clocks (low-mass, low power-consumption, frequency 



O ' stability comparable to that of ground-based clocks) can enable interplanetary spacecraft radio 

Q ■ science experiments at unprecedented Doppler sensitivities. The addition of an on-board digital 

receiver would allow the up- and down-link Doppler frequencies to be measured separately. Such 
separate, high-quality measurements allow optimal data combinations that suppress the currently- 
leading noise sources: phase scintillation noise from the Earth's atmosphere and Doppler noise 



caused by mechanical vibrations of the ground antenna. Here we provide a general expression 
for the optimal combination of ground and on-board Doppler data and compute the sensitivity 
such a system would have to low-frequency gravitational waves (GWs). Assuming a plasma scin- 

00 

■ tillation noise calibration comparable to that already demonstrated with the multi-link CASSINI 

(N 

. radio system, the space-clock/digital-receiver instrumentation enhancements would give GW strain 

OO sensitivity of 2.0 x 10 for randomly polarized, monochromatic GW signals over a two-decade 

o 

(~ 0.0001 — 0.01 Hz) region of the low-frequency band. This is about an order of magnitude better 
^ ■ than currently achieved with traditional two-way coherent Doppler experiments. The utility of 

H : 

. . , optimally combining simultaneous up- and down-link observations is not limited to GW searches. 

The Doppler tracking technique discussed here could be performed at minimal incremental cost 
to also improve other radio science experiments (i.e. tests of relativistic gravity, planetary and 
satellite gravity field measurements, atmospheric and ring occultations) on future interplanetary 
missions. 
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I. INTRODUCTION 



Measurements of the relative velocity between the Earth and an interplanetary spacecraft, 
by ,„ea,. of coherent microwave tracking, have allowed stud.es of solar system bodies g 



tests of relativistic gravity , searches for low-frequency gravitational radiation \i, U 



and other science objectives [1]. In the frequency band (10^^ — 10~^) Hz, typical deep space 
tracks are limited by phase scintillation caused by random refractivity variations in the solar 
wind, and the ionosphere [?]. The most sensitive deep-space Doppler observations to date, 
however, calibrate and largely remove these noises [3|, Il0| and are then limited by 

antenna mechanical noise (unmodeled motion of the phase center of the ground antenna) 
and residual post-calibration tropospheric scintillation (i.e. Doppler fluctuations caused by 
refractive index fluctuations in the Earth's atmosphere) The most sensitive obser- 

vations hit the limit identified by these noise sources with an Allan standard deviation of 
about 3 X 10~^^ for integration times of a few thousand seconds. Improved sensitivity would 
benefit the science disciplines listed above, but antenna mechanical noise, in particular, has 
seemed irreducible at reasonable cost since it would require a large, moving, steel structure 
much more rigid than that of the current ground tracking stations. 

Several ideas have been proposed to reduce the antenna mechanical noise or the tro- 



pospheric noise, or both 



12| . Those ideas do not involve modifications to the spacecraft. 



but rather use simultaneous tracking on the ground with appropriate linear combinations of 
those data to synthesize an observable with the ground/tropospheric noises of the better of 
the two receiving stations. If, however, additional microwave instrumentation on th e sp ace- 



craft is considered, such as a space- 
a space-quahfied digital receiver [15, 



3orne, high-stability frequency standard 
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IJ] and 



161 ]. then an alternative method for noise suppression 



is possible. This involves properly combining the two one-way (spacecraft to Earth and 
simultaneous Earth to spacecraft) Doppler data taken onboard and on the ground in such 
a way to maximize the signal-to- noise ratio of the observed physical observable.^ Multiple 
Doppler observations with different noises, or different transfer functions to the same noises. 



^ It should be emphasized that the tropospheric, antenna mechanical and ionospheric noise suppression 
can also be accomplished by combining the two-way coherent Doppler data with the one-way Doppler 
measurement performed at the ground [6|,ll9|. We have however analyzed the configuration with the two 
one-way measurements for reasons of symmetry and simplicity. 
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are clearly useful in identifying noise sources and minimizing (and in some cases canceling) 
their effects on the final observable. In particular, the use of multiple radio links (some of 
which driven by a high-quality space-borne fre quen cy standard) was pioneered by R. Vessot 
in the Gravity Probe A sub-orbital experiment jl7|]. Because of mass and power constraints, 
however, no very high quality frequency standards have yet been flown on deep space probes. 

Recent advances in clock technology indicate that a new era of space-qualified, highly 
stable frequency standards has started [isl, which will result into significantly improved 
Doppler radio science experiments. A summary of this paper is given below. 

In Section [III we present a brief overview of the theory underlying Doppler tracking ex- 
periments relying only on the two one-way Doppler data measured onboard and on the 
ground. Although this experimental configuration has been discussed in previous publica- 
tions [19 1, it has been shown only relatively recently 0, Q] how to fully take advantage 
of it for improving the precision of Doppler tracking radio science experiments. A brief de- 
scription of an onboard atomic clock and a digital receiver is given in the Appendix, where 
an account of all the one-sided power spectral densities of the noises affecting the Doppler 
data is also presented. 

The main advantage of spacecraft Doppler tracking experiments relying on the two one- 
way Doppler data, over those based on two-way coherent measurements, is in their ability of 
exactly canceling the frequency fluctuations due to the Earth atmosphere and ionosphere, 
and the mechanical vibrations of the ground antenna, presently the main noise sources 

n 

of Doppler tracking experiments [6|. This is possible because there exists a unique linear 
combination of the properly time-shifted one-way measurements that does exactly that Q]. 
Depending on the specific radio science experiment performed with this technique, it is actu- 
ally possible to combine optimally the two one-way Doppler measurements to maximize the 
signal-to-noise ratio (SNR) of the experiment. After deriving the expression of the optimal 
SNR in Section UTTl we apply it to searches for gravitational waves. Under the assumption of 
calibrating the frequency fluctuations induced by the interplanetary plasma, we find that a 
Doppler broad-band sensitivity of 2.0 x 10^^^ to randomly polarized monochromatic signals 
uniformly distributed over the sky can be achieved. This is about one order of magnitude 
better than that obtainable with two-way coherent Doppler tracking experiments. Narrow- 
band searches at frequencies where the transfer function of the onboard clock reaches sharp 
nulls (i.e. the "xylophone" frequencies) Gj further enhance the strain sensitivity of these 
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experiments to about 7.0 x 10~^^. 

In Section IIVI we finally present our comments and conclusions, and emphasize that the 
Doppler tracking technique discussed in this article can be performed at minimal additional 
cost by forthcoming interplanetary missions. 



II. THE ONE-WAY DOPPLER TRACKING OBSERVABLES 



In Doppler tracking experiments aimed at detecting low-frequency (milliHertz) gravita- 
tional radiation, a distant interplanetary spacecraft is monitored from Earth through a radio 
link, with the Earth and the spacecraft acting as free-falling test particles ^. A radio signal 
of nominal frequency z/q is transmitted to the spacecraft, and coherently transponded back 
to Earth where the received signal is compared to a signal referenced to a highly stable 
clock (typically a hydrogen-maser) . Relative frequency changes Azz/z/q, as functions of time, 
are measured. When a gravitational wave crossing the solar system propagates through the 
radio link, it causes small perturbations in Azz/z/q, which are replicated three times in the 
Doppler data with maximum spacing given by the two-way light propagation time between 
the Earth and the spacecraft ^. 

An alternative way of performing Doppler tracking searches for gravitational radiation 



was suggested in 



191]. By adding a highly-stable clock and a digital receiver to 
the spacecraft radio instrumentation (Figure [1]), two one-way Doppler time series can be 
recorded simultaneously at the ground station and on board the spacecraft. 

If we introduce a set of Cartesian orthogonal coordinates {X, Y, Z) in which the wave 
is propagating along the Z-axis and (X, Y) are two orthogonal axes in the plane of the 
wave (see Figure [2]), then the two one-way relative frequency fluctuations at time t can be 
written in the following form after first-order Doppler and other systematic Doppler effects 
are modeled out from the data ^|6| 



^ Spacecraft Doppler GW searches piggy-back on interplanetary probes used primarily for other (e.g., solar 
system) science goals. Doppler tracking is the current generation GW detector in the low- frequency 
band. A much more sensitive, dedicated GW mission, LISA, is currently in the design and technology 
development stage and could launch sometimes in the next decade. 

The one-way Doppler data measured onboard can be digitally recorded, time tagged, and telemetered 
back to Earth in real time or at a later time during the mission. 
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FIG. 1: Block diagram of the radio hardware at the ground antenna of the NASA Deep Space 
Network (DSN) and on board the spacecraft (S/C), which ahows the acquisition and recording of 
the two Doppler data E{t), S{t). A description of each individual block in this diagram is provided 
in the Appendix. 



V ^0 A 



[h{t-{l + fi)L) - h{t)] + Cs{t-L) - CE{t) 



+ T{t) + B{t-L) + As{t-L) + ELE{t) + P^;(i) , (1) 



^""^^^^ -S{t) = ^ [h{t-L) - h{t-iiL)] + CE{t-L) - Cs{t) 



1^0 



+ T{t -L) + B{t) + AE{t -L) + ELs{t) + Ps{t) , (2) 



where h{t) is equal to 



h(t) = h+{t) cos(2(/)) + /ix {t) sin(2(/)) . 
6 



(3) 




FIG. 2: A radio signal of nominal frequency is transmitted to a spacecraft and simultaneously 
another radio signal from the spacecraft and referenced to the onboard clock is transmitted to the 
ground. The gravitational wave train propagates along the Z direction, and the cosine of the angle 
between its direction of propagation and the radio beam is denoted by /i. See text for a complete 
description. 

Here h+(t), hx{t) are the wave's two amplitudes with respect to the {X,Y) axis, (6', 0) 
are the polar angles describing the location of the spacecraft with respect to the {X, Y, Z) 
coordinates, fi is equal to cos^, and L is the distance to the spacecraft (units in which the 
speed of light c = 1) 

In Eqs. ([H [2]) we have assumed the Earth and the on board clocks to be perfectly synchro- 
nized. Although this condition is impossible to achieve in practice, it has been previously 
shown by one of us ^ that the accuracy required for successfully implementing a noise can- 
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cellation scheme similar to the one discussed in this paper requires a clock synchronization 
accuracy of about 0.5 s, which is easy to achieve. 

In Eqs. ( TTPj) we have denoted by CE(t), Cs(t) the random processes associated with the 
frequency fluctuations of the clock on Earth and onboard respectively, B{t) the joint effect 
of the noise from buffeting of the probe by non gravitational forces and from the antenna 
of the spacecraft, T(t) the joint frequency fluctuations due to the troposphere, ionosphere 
and ground antenna, Asit) the noise of the radio transmitter on the ground, As{t) the noise 
of the radio transmitter on board, ELE{t), ELs{t), the noise from the electronics on the 
ground and onboard respectively, and PE(t), Ps(t) the fluctuations on the two links due 
to the interplanetary plasma. Since the frequency fluctuations induced by the plasma are, 
to first order, inversely proportional to the square of the radio frequency, by using high 
frequency radio signals or by monitoring two different radio frequencies transmitted to and 
from the spacecraft, this noise source can be suppressed to very low levels or entirely removed 

n 

from the data respectively [9']. In what follows we will assume dual frequency to be used, 
and disregard the noise effects of the plasma fluctuations in our analysis. 

From Eqs. ( |T|2|) we deduce that gravitational wave pulses of duration longer than the one- 
light-time L give a Doppler response that, to first order, tends to zero. The tracking system 
essentially acts band device, in which the low- frequency limit /; is roughly equal to 

Hz, and the high-frequency limit fn is set by the thermal noise in the receiver. Since 
the clocks and some electronic components are most stable at integration times around 1000 
seconds, Doppler tracking experiments are performed when the distance to the spacecraft is 
of the order of a few astronomical units. This sets the value of /; for a typical experiment 
to about 10~^ Hz, while the thermal noise gives an fn of about 10~^ Hz. 

It is important to note the characteristic time signatures of the clock noises, CE{t) and 
Cs(t), of the probe antenna and spacecraft buffeting noise B(t), of the troposphere, iono- 
sphere, and ground antenna noise T(t), and the transmitters AE{t), As{t). The time signa- 
ture of the two clock noises, for instance, can be understood by observing that the frequency 
of the signal received at the ground station at time t contains fluctuations from the onboard 
clock that were transmitted L seconds earlier and the noise from the ground clock enters 
with a negative sign at time t due to the heterodyne nature of the Doppler measurement. 
The time signature of the noises T, B(t), AEit), and As(t) in Eq. ([1121) can be understood 
through similar considerations. 
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Since the major noise source affecting the two one-way measurements is represented 
by the ffuctuations induced by the Earth troposphere and the mechanical vibrations of the 
ground station, it has been emphasized gI that there exists a combination of the two Doppler 
data that cancels these noises. It is easy to see from inspection of Eqs. f[T|2]) that such a 
combination is equal to 

x{t) = S{t) - E{t - L) . (4) 
After substituting into Eq. (jl]) the expressions for -E(t), S{t) given in Eqs. (fT|2l) we get 

x{t) = h{t -L)- h{t - liV) - h{t -2L- fiL) 

+ 2CE{t -L)- [Cs{t) + Cs{t - 2L)] + [B{t) - B{t - 2L)] 

+ AE{t -L)- As{t - 2L) + ELs{t) - ELE{t - L) . (5) 

From Eq. ([5]) we may notice that the spacecraft buffeting noise, B, does not cancel 
exactly and it gets suppressed by its transfer function to the x combination at frequencies 
smaller than the inverse of the round-trip-light-time. 



III. GRAVITATIONAL WAVE SENSITIVITIES 



This section describes the derivation of the sensitivity, which is defined on average over 
the sky, to be equal to the strength of a sinusoidal gravitational wave required to achieve a 
signal-to-noise ratio of 1 in a forty-day integration time. Note that the sensitivity is therefore 
a function of Fourier frequency, /. We have chosen the integration time to be equal to forty 
days since this was the tracking time of the CASSINI gravitational wave experiments |5|. 
Sensitivity is essentially the noise-to-signal ratio and it will be computed for both the new 
data combination x as well as for the traditional two-way coherent tracking measurement, 
y, for comparison reasons. For convenience we provide below the expression of the two-way 
Doppler response, y{t), which will be used later on for estimating its sensitivity 



yit) = -^^^ Ht) - fihit-{l+ + ^i±i^ h{t - 2L) 
+ CE{t-2L) - CEit) + 2B{t-L) + T{t-2L) + T{t) 
+ AE{t-2L) + Asit-L) + TR{t-L) + ELE,{t) , (6) 
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where TR is the random process associated with the relative frequency fluctuations due to 
the onboard microwave transponder. For more details we refer the reader to 21|. 



A. Signal Averaged Power 

The averaged signal power in the combination x, estimated at an arbitrary Fourier fre- 
quency /, is computed by (i) taking the Fourier transform of the signal entering into the 
combination x and its modulus-squared, and (ii) by integrating the resulting expression 
over an ensemble of sinusoidal signals uniformly distributed over the celestial sphere. Such 
a calculation is long but straightforward, and the resulting expression, S'^;^, is equal to 



+ - cos (ZTT/L) 



(7) 



3 3 ^ ■> ^ 2(7r/L)2 

where Sh is the gravitational wave signal one-sided power spectral density. 

The calculation of the averaged signal power of the observable y can similarly be carried 
through, resulting into the following expression 

Sn'fL' + 2 sin(47r/L) - |7r/L (3 + An'fL^) cosjAn fL) - QnfL 

Figure [3] shows the two "transfer functions" of the signal one-sided power spectral den- 
sities, i.e. = Sx^/Sh and qy = Sy^/Sh- The transfer function is slightly larger than 
qy in the region [5 x 10^^ — 1] Hz, indicating that a constructing interference of the signal 
with itself is taking place in this part of frequency band. On the other hand, in the low-part 
of the frequency band the combination x shows a coupling to gravitational radiation that 
is weaker than that of y. This is to be expected, as x is the difference of the two one-way 
measurements, which in the "long-wavelength" limit become equal to each other. 



B. Noise Spectra 

To compute the sensitivity of the combination x, and compare it against that of the 
two-way Doppler measurement, y, we need the one-sided power spectral densities of the 
main noise sources and their transfer functions to the observables x and y. If we assume 
all the noise sources to be uncorrelated, from equations ([5], [6]) we can derive the following 
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Signal Transfer Functions 




FIG. 3: Averaged power transfer functions of the Doppler responses x and y to an ensemble of 
sinusoidal signals randomly polarized and uniformly distributed over the celestial sphere. The 
x-transfer function shows constructing interference at frequencies that are integer multiples of the 
inverse of the round-trip- light-time, 2L, taken here to be equal to 5500 seconds. The coupling of 
the X data combination to such a stochastic ensemble of gravitational radiation is slightly stronger 
than that of the two-way y response at frequencies larger than 5.0 x 10~^ Hz. At lower frequencies 
the transfer function of the x combination decays more rapidly than that of y as a consequence of 
being the difference of the two one-way measurements. 

expressions for the two noise spectra Sx„ and Sy^ 

= 4Sc, + ^Scs cos2(27r/L) + 4Sb sm\27T fL) + Sa, + + Sel, + Sels , (9) 
Sy„ = ASc^ sin2(27r/L) + ASt cos\2n fL) + ASb + Sa^ + Sas + Sel^ + Sels + StrHO) 

where the meaning of the various terms appearing into Eqs. fpiTU]) is self-explanatory. We 
provide in the Appendix the expressions for the various noise spectra corresponding to a 
gravitational wave search performed with a spacecraft out to a distance of 5.5 AU from Earth 
(as was during the first gravitational wave experiment with the CASSINI spacecraft), and 
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FIG. 4: Sensitivities of the x (solid-line) and y (dash-line) Doppler responses to a randomly polar- 
ized and uniformly distributed stochastic ensemble of sinusoidal gravitational wave signals. The 
sensitivity is expressed as a function of the frequency and represents the equivalent sinusoidal strain 
required to produce a signal-to-noise ration of 1. The spacecraft is assumed to be out to a distance 
of 5.5 AU, and eighty percent tropospheric noise calibration is applied to the two-way Doppler data 
y (dash-line). The two sensitivity curves reflect the noise spectral levels and shapes (given in the 
Appendix), their transfer functions to the observables x and y (Eqs. [9l llOp . and the gravitational 
wave transfer functions shown in Fig. [3]) 

equipped with an onboard microwave instrumentation similar to the one flown on CASSINI. 

The gravitational wave sensitivity is the wave amplitude required to achieve a signal- 
to- noise ratio of 1, and it can be computed as a function of Fourier frequency using the 
following expression [Sj 

where z means either x or y, and B is the frequency bandwidth corresponding to a 40 days 
integration time, the duration of the gravitational wave experiments performed with the 
CASSINI spacecraft fl. 
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The main difference between the two observables x and y is of course in the absence in the 
X combination of the joint disturbances from the troposphere and the mechanical vibration 
of the ground antenna (the random process denoted T in Eq.([2])) and the presence (in x) of 
the spacecraft clock noise process. As onboard and ground microwave instrumentation have 
in recent years reached unprecedented frequency stabilities, T has indeed become the major 
sensitivity limitation of two-way Doppler tracking searches for gravitational radiation |5|. 
Since the effects from the troposphere can be mitigated by relying on simultaneous mea- 
surements performed by a radiometer located in the proximity of the tracking station, our 
sensitivity analysis reflects the assumption of being able to calibrate out eighty percent of the 
tropospheric effects from the y observable (as demonstrated by the CASSINI experiments). 
Figure H] shows the estimated sensitivity of the observable x (continuous-line) formed out of 
the two one-way measurements, and compares it against that of the two-way measurement 
in which eighty percent of the tropospheric effects are calibrated out (dash- line). The x 
combination displays the best sensitivity in the frequency band 10~^ — 10^^ Hz, which is 
of most interest to gravitational wave search experiments. At higher frequencies, between 
about 10"^ — 3.0 X 10^^ Hz, effects related to the locking of the atomic clock to its local 
oscillator introduce a small degradation in the x sensitivity over that of the y data. Also, 
at frequencies lower than 3.0 x 10~^ Hz, the x combination shows a sensitivity worse than 
that of y due to the cancellation of the signal in this low-frequency region. 

As the two one-way Doppler measurements can be combined to synthesize the two-way 
measurement y {gI, Q , one could argue that at these frequencies one could of course rely 
on the synthesized y data to take advantage, if needed, of its better sensitivity in these two 
regions of the accessible frequency band. This observation suggests that it must be possible 
to identify a combination of the two one-way Doppler data that maximizes the sensitivity 
to gravitational waves in the entire band of interest. 



C. Optimal Sensitivity 

In order to derive the combination of the two one-way Doppler data that achieves optimal 
sensitivity, let us consider the following linear combination ?](/) of the Fourier transforms 
of E{t) and Sit) 

r/(/) ^ ai(/. A) EU) + a^U. A) 5(/) , (12) 
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where the {ai{f, A)}j=i^2 are arbitrary complex functions of the Fourier frequency /, and of 
a vector A containing parameters characterizing the signal and the noises affecting the two 
Doppler data. For a given choice of the two functions {aj}j=i_2, rj gives a specific Doppler 

ying, for a given signal, the two 



26( 1 ■ SNR'i, of the combination rj 



data combination, and our goal is therefore that of identi 
functions {cii}i=i,2 that maximize the signal-to- noise ratio 

SNg}, = + df . (13) 

Jfl {\ai En + a2 SnP) 

In equation f|T3|) the subscripts s and n refer to the signal and the noise parts of {E, S) respec- 
tively, the angle brackets represent noise ensemble averages, and the interval of integration 
ifh fu) corresponds to the accessible frequency band. 

The SNR?^ can be regarded as a functional over the space of the two complex functions 
{cii}i=i,2, and their expressions that maximize it can of course be derived by solving the 
associated set of Euler-Lagrange equations. The derivation of the expression of the optimal 
SNR, SNR"^^^^, is long but straightforward, and it is equal to (see [27| for details) 

SNR' = rzp(C-^),,zf rf/. (14) 



In Eq. ( fT4ll the convention of sum over repeated indices is assumed, z^*-* is the vector of 
the signals, [Es^ Sg), and C is the hermitian, non-singular, correlation matrix of the vector 
random process = (^E^, Sn) 

(CU ^ (z(")zi">) . (15) 



Eq. f[T4|) can now be used for estimating the sensitivity to an ensemble of sinusoidal 
gravitational wave signals randomly polarized and uniformly distributed over the celestial 
sphere. Figure [5] shows the estimated optimal sensitivity (solid-line) obtained by relying on 
the two one-way measurements, and again that of a two-way coherent tracking experiment 
(dash-line), in which all the parameters characterizing the experiment are as in Figure (jll). 
Note how the sensitivity of the optimal combination is now consistently below that of the 
two-way measurement, and it coincides with that of the x combination in most of the 
accessible frequency band. 
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Sensitivities 




FIG. 5: Sensitivity curves of the optimal combination (solid-line) and the two-way Doppler tracking 
data (dash- line). The noise parameters are equal to those used in figured! Note how the sensitivity 
of the optimal combination is now consistently below that of the two-way measurement. 

IV. CONCLUSIONS 

We have discussed a method for significantly enhancing the sensitivity of Doppler tracking 
experiments aimed at the detection of gravitational waves. The main result of our analysis 
shows that by adding an atomic frequency standard and a digital receiver on board the 
spacecraft we can achieve a broad-band sensitivity of 2.0 x 10^^'' in the milliHertz band. 
This sensitivity figure is obtained by completely removing the frequency fluctuations due to 
the interplanetary plasma. Our method relies on a properly chosen linear combination of the 
one-way Doppler data recorded on board with the data measured on the ground. It allows 
us to optimally suppress the frequency fluctuations due to the troposphere, ionosphere, and 
antenna mechanical and, for a spacecraft that is tracked for 40 days out to 5.5 AU, to reach a 
sensitivity that is about one order of magnitude better than that achievable by a state-of-the 
art two-way Doppler tracking search. 

The expression of the optimal combination of the two one-way Doppler data can be used in 
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all the classic tests of relativistic theory of gravity in which one-way and two-way spacecraft 
Doppler measurements are used as primary data sets. We will analyze the implications of 
the sensitivity improvements that this technique will provide for direct measurements of the 
gravitational red-shift, the second-order relativistic Doppler effect predicted by the theory 
of special relativity, searches for possible anisotropy in the velocity of light, measurements 
of the parameterized post-Newtonian parameters, and measurements of the deflection and 
time delay by the Sun in radio signals. This research is in progress, and will be the subject 
of a forthcoming investigation. 
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Appendix 

A. Noise sources and their spectra 

This appendix provides a general description of the radio hardware needed for imple- 
menting the technique described in the main body of this paper, the corresponding one-sided 
power spectral densities of the frequency fluctuations introduced by these subsystems into 
the observables E{t) and S(t), and discusses the frequency fluctuations due to the Earth 
atmosphere. For a more comprehensive analysis on the radio hardware the reader is referred 
to 



aospnere. i^or a more comprenensive analysis on tne radio nardware tne readt 

Q, y, Q, Q, Q, Q, while a review on the propagation noises is given in fl 



The ground master clock and the frequency & timing distribution represent the overall 
contribution of the reference clock itself and the cabling system that takes the signal gener- 
ated by the master clock to the antenna. This can be located several kilometers away from 
the site of the clock, implying that the need of a highly-stable cabling system is required. 
It has been shown at JPL that optical-fiber cables would not degrade significantly the fre- 
quency stability of the signal generated by the master clock. The corresponding one-sided 
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power spectral 



is equal to 



22| 



density of the frequency fluctuations, introduced by these two noise sources, 



Sc^if) = 6.2 X [lO'^V + 10-^^ f-^ + 10"^°] + 1.3 X lO-^V^ Hz"^ . (16) 



The Ground and onboard Transmitters include all the frequency multipliers that are 
needed to generate the desired frequency of the transmitted radio signal, starting from the 
frequency provided by the clocks. It also accounts for the radio amplifier, and the extra 
phase delay changes occurring between the amplifier and the feed cone of the antennas. The 
noise due to the amplifiers is the dominant one, and it has been characterized in j^. The 
one-sided power spectral densities of the frequency fluctuations are given by 

SaAI) + SAsif) = 2.3 X 10-28 + 4.0 X lO-^V Hz"^ , (17) 

The noises introduced by the Receiver at the ground station can be modeled as white 
phase fluctuations. The contribution to the overall noise budget from the receiver chain on 
the ground can be repartitioned into thermal noise (finiteness of the signal-to-noise ratio) and 
fluctuations introduced into the signal as it propagates through the cables and waveguides 
running from the feed of the antenna to the actual receiver. The effects of the latter noise 
source is nowadays minimized with the use of beam waveguide (BWG) antennas. These new 
antennas have become operational in the year 2004 at the NASA Deep Space Network three 
sites: in North America (Goldstone, California), Europe (Madrid, Spain), and Australia 
(Canberra). 

Under the assumption of relying on a 34 meter diameter beam waveguide antenna for 
receiving a coherent Ka-Band (32 GHz) signal transmitted by a spacecraft out a distance of 
5.5 AU, a ground system noise temperature of about 70 degrees Kelvin, an onboard Ka-Band 
amplifier of 10 W, and a spacecraft High Gain Antenna (HGA) with a diameter of about 4 
meters, we find the following one-sided power spectral density of the frequency fluctuations 
at Ka-Band 

5'i^L,(/) = 6.3 X 10-2 V'Hz-^ (18) 

The buffeting of the spacecraft will introduce unwanted frequency fluctuations in the one-way 

I I 

Doppler observable. Estimates of its magnitude have been given in ^23|], and the one-sided 
power spectral density of the frequency fluctuations is given by the following expression 

5b(/) = 5.0 X 10-^2 f-S ^ ^ jj^-l _ (^g) 
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The noise introduced into the Doppler observable y by the Earth Atmosphere and ionosphere, 
and by the scintillation of the interplanetary plasma, have been studied extensively in the 
literature 



24|. 



The scintillation introduced into the Doppler observables by the Atmosphere are indepen- 
dent of the microwave frequency at which the spacecraft is tracked. Since gravitational wave 
searches are performed in a band whose upper frequency cutoff is smaller than 1 Hz (ther- 
mal noise at higher frequencies becomes unacceptably large), the one-sided power spectral 
density of the noise due to the atmosphere can be written as follows {25! 



ST{f) = 2.8 X 10-2^ /-^/^ Hz"^ 10"^ < / < 10-2 Hz 

= 2.2 X 10"^° f-^ Rz-^ 10^2 < / < 1 Hz . (20) 

The first term in the equation above accounts for the remaining effect of the atmosphere after 
eighty percent calibration is applied to the data with the use of a water vapor radiometer , 
while the second term accounts for the effect of aperture averaging, which causes a reduction 
in delay fluctuations on time scales less than the antenna wind speed crossing time (1 to 10 
seconds) {25^. 

The Transponder entering into the y measurement is responsible for keeping the phase 
coherence between the incoming and outgoing radio signals on the spacecraft. Its perfor- 
mance depends on the accuracy of tracking of the up-link signal b y th e phase locked loop, 
and the noise floor and non-linearities of its electronic components 3]. Frequency stability 
measurements of the Ka-Band (32 GHz) transponder flown onboard the CASSINI mission 
have resulted in the following one-sided power spectral density of the relative frequency 
fluctuations 

Smif) = 1.6 X 10-26 / Hz-i . (21) 

The onboard clock provides the frequency and timing reference for the onboard radio in- 
strumentation, and identifies the stability of the microwave signal transmitted to the ground. 
The space-qualified clock presently under realization at the Jet Propulsion Laboratory relies 
on a combined "interplay" between a local quartz oscillator and a trapped Hg ions clock. The 
frequency of the hyperfine transition made by the Hg ions is used for monitoring and cor- 
recting the frequency of the quartz oscillator. This steering process takes place on a typical 
time scale of about 10 seconds, making then possible over longer time scale to significantly 

18 



improve the stability of the resulting combined instrument. A frequency stability compara- 
ble to that of the Deep Space Network ground clocks has already been demonstrated with 
a prototype, showing an Allan standard deviation of a few parts in 10~^^ at an integration 



time of a few thousand seconds 



13 



14| . The corresponding one-sided power spectral density 



of the relative frequency fluctuations of such a clock is given by the following expression 



10"^ < / < 2.0 X 10-^ Hz 
2.0 X 10"^ < / < 2.0 X 10"^ Hz 
2.0 X 10"^ < / < 1 Hz (22) 

The onboard digital receiver used in this method measures the amplitude and phase of 
the uplink signal and telemeters that information to the ground. 

The frequency fluctuations of the receiver chain on board the spacecraft is estimated to be 
entirely due to thermal noise because of a simpler cabling system, and its performance being 
essentially identical to that of the ground receiver. By assuming again a 34 meter diameter 
beam waveguide antenna transmitting with an 800 W Ka-Band (32 GHz) amplifier to a 
spacecraft out to a distance of 5.5 AU, equipped with a 4 meter diameter (HGA), and 
a system noise temperature of 400 Kelvin, we find the following one-sided power spectral 
density of the frequency fluctuations of the onboard electronics noise at Ka-Band [sl 

SelsU) = 7.2 X 10-^' fRz-\ (23) 



Scs - 



= 5.0 X 10"^^ Hz"^ 
= 2.5 X lO^^V Hz~^ 
= 10^26^-1 Hz-i 
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